Factorial moment analyses in difFractive lepton-nucleon scattering 



Zhang Yang 

Institut fiir theoretische Physik, FU Berlin, Arnimallee I4, 14195 Berlin, Germany 



Abstract 

It is pointed out that "the colorless objects" in diffractive lepton-nucleon scattering 
in the small-xs region can be probed by measuring the scaled factorial moments of 
final-state-hadrons and the dependence of their scaling behavior upon the diffractive 
kinematic variables. The Monte Carlo implementation of RAPGAP and JETSET are 
discussed as illustrative examples. The results of these model-calculations show in 
particular that inclusion of the contributions from the gamma gluon fusion processes 
can considerably enlarge the power of the scaled factorial moments. The possibility 
for probing the anomalous scaling behaviors of probability moments of the transverse 
energies in HERA calorimeter environment is also discussed. 







Recent HERA-experiments on deep-inelastic electron-proton (ep) scattering (DIS) in the 
low xb kinematic range have clearly shown the existence of a distinct class of events. These 
events are characterized by the fact that there is no hadronic energy flow in a considerably 
large interval of pseudo-rapidity t] adjacent to the proton beam direction. Our present 
understanding of DIS could be inadequate at low-x^ because additions to the leading order 
QCD-based partonic picture are likely to be substantial. A natural interpretation of these so 
called "large rapidity gap" events is based on the hypothesis that the deep-inelastic scattering 
process involves the interaction of the virtual boson probe with a colorless component of 
the proton. Hence there is no chromodynamic radiation in final state immediately adjacent 
to the direction of the scattered proton or any proton remnant. What is this colorless 
component originating from the nucleon? The large rapidity gap events discovered in deep 
inelastic scattering at HERA0 are usually interpreted in terms of pomeron-exchanged mo dei. 
Although this seems to work reasonably well phenomenologically, there is yet no satisfactory 
understanding of the pomeron structure and its interactions mechanism. In this respect, it 
is helpful to probe the properties of such objects in the non-traditional aspects, in addition 
to the traditional measures such as rapidity gap distribution, structure function and the 
averaged cross section and so on0. In the present note we wish to point out that useful 
information about the exchanged colorless object in diffractive lepton-nucleon scattering 
process can be extracted by studying scaling behavior and fractality (intermittency) of the 
final state of the colorless component in proton excited by the virtual boson probe. 

The manifestation of fractality (intermittency) in high energy multiparticle production 
process is the anomalous scaling!'! 

/ Ax \ '^'^ 

Fg{6x) = Fg{Ax) ( — ] as M oo, 6x (1) 

of g-order factorial moments (FM's) Fg, defined as 

where x is some phase space variable, e.g. (pseudo-)rapidity, the scale 6x = Ax/M is the 
bin width for a M-partition of the region Ax in consideration, is the multiplicity in the 



mth bin. Since the factorial moments can remove the statistic noise around the probabihty 
and associated directly with the scaled moments of probabilities!, the scaling exponent (l)q 
in Eq.( |1|), called intermittency index, characterizing the strength of dynamical fluctuation, 
is connected^ with the anomalous fractal dimension dq of rank q of the spatial-temporal 
evolution of high energy collisions, dq = (pq/{q — 1). A possible cause leading to the power- 
law of FM's of final particles in high-energy collisions is that the emitting source of final 
hadrons is self-similar fractali. Furthermore, DIS experiments and the empirical analysesBB 
show that the gluon-density in the nucleon in the low-x^ kinematical region is much higher 
than those for quarks/ ant iquarks, and it is increasing with decreasing x^B. In this soft gluon 
system where the gluons interact with each other in complicated dynamic processes, it has 
been proposedll that the gluons may 'self-organize' into the clusters in the dissipative gluon 
system by self-organized criticality (SOC)i, and the colorless component in proton can be 
regarded as color singlet gluon cluster^^l. The spatial-temporal structure of the SOC-cluster 
(or BTW-cluster) is self-similar fractali. So it is feasible to study the fractal structure of the 
colorless component of the proton by measuring the anomalous scaling behaviors of factorial 
moments of the final state particles originating from the scattering of virtual photon and 
colorless object in the diffr active lepton-nucleon scattering. 

The xs-dependence of scaling behaviors of FM's, and in particular that of intermittency 
index (pq{Q^,XB',xp) for fixed (and xp), plays a distinguished role in such studies. This 
is because, when the virtual photon 7* originating from the incident lepton is absorbed 
by the nucleon, its energy- momentum q = {q^,q) is in fact absorbed by a virtual colorless 
component of the nucleon in the lepton-nucleon diffractive scattering. In a fast moving 
reference frame, for example the lepton-nucleon center-of-mass frame, where the nucleon's 
momentum P is large in high-energy collisions, the time interval Tint in which the absorption 
process takes place (it is known as the lepton-nucleon interaction/collision time) can be 
estimated by making use of the uncertainty principle. In fact, by calculating 1/q^ in this 
reference frame we obtainS: 
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This means, for given P and = —q"^, the interaction time Tint is directly proportional 
to in the low-x^ {xb ^ 1) kinematic range. In other words, is a measure of the 
time-interval in which the absorption of 7* by the space-like virtual colorless object takes 
place. Hence, by studying the x^-dependence of the intermittency index x^; xp), 

we are not merely probing the anomalous scaling behaviors of the collision process between 
7* and the colorless object which we hereafter call Cg. Since this hadronization process of 
the virtual colorless object Cg is initiated by the interaction with 7*, we are also examing 
whether/how the hadronization process changes with the interaction time Tint- This question 
is of considerable interest, because a virtual photon 7* can (logically) only be absorbed 
by virtual systems (cg's) whose lifetimes (t*'s) are longer than interactive time Tint (i.e. 
Tint ^ Tc). That is, the average lifetime (tc) of the Cq's, which can absorb a 7* associated 
with interaction-time Tint, is a function of Tint- Hence, from the x^-dependence of the scaling 
behavior of FM's, in particular from that of the corresponding ^^(xb, Q^; xp)'s, we can in 
principle find out whether/how the dynamic fluctuation and the fractal structure in the 
spatial-temporal evolution originating from the exchanged colorless object Cq depends on its 
average lifetime (tc) of the Cq's. 

The Q^-dependence of the scaling behavior of the FM's is also of considerable interest. 
This is because, in photon-proton scattering experiments, not only those with real (Q^ = 
0) photons but also those with space-like (Q^ > 0) photons where is not too large 
(< lGeVVc^ say) have very much in common with hadron-hadron collisions. Having in 
mind that the index of intermittency for hadron-hadron scattering is smaller than that for 
electron-positron annihilation processes^, we are led to the following questions: Do we expect 
to see a stronger Q^-dependence when we increase from zero to 10 or 100 GeV'^/c^, say? 
Is this also a way to see whether space- like photons at large "behave like hadrons" in 
such interactions? 

While waiting for data to perform the above-mentioned analysis, let us now consider the 



following two plienomenological models as illustrative examples: 

(A). If the colorless object (cq) is a quark-antiquark pair (formed by interacting gluons) 
which exists in the time-interval when the virtual photon 7* is absorbed by the object, 
we shall see the following: Especially when is sufficiently large, the incoming 7* (the 
transverse dimension is expected to be proportional to l/Q^) will hit the quark (g) or the 
antiquark (g) and make them fly apart symmetrically with respect to the center of mass 
the 7*gg system — similar to the gg-pair produced in e"'"e~-collisions (with respect to the 
center of mass the gg system). That is, in this case, the flnal-state-hadrons of an event are 
fragmentation-products of the quark and/or the antiquark, and hence they are expected to 
show characteristic features similar as those observed in the reaction e^e" hadrons at 
the same c.m.s. energy. It is interesting to see that the very recent inclusive measurements 
performed at heraS in diffractive electron-proton scattering show the following: Both 
the scaled longitudinal momentum [xp) distribution and the transverse momentum (p^) 
distribution are strikingly symmetric with respect to the center-of-mass of the photon and 
the struck colorless object; and the general features of these distributions are very much the 
same as those observed in electron-position collision processes. These facts strongly suggest 
that a more detailed comparison between these two collision processes would be useful. 

For this purpose, we made use of the Monte Carlo (MC) program JETSEtO. We 
generated 50,000 MC events, and calculated the second normalized factorial moment in 3- 
dimensional (77, 0) phase space at the given gg-c.m.s energy i/i, where the pseudorapidity 
7], transverse momentum p_L and the azimuthal angle are deflned with respect to the 
sphericity axis of the event. The usual cumulative variables X translated from x = [r], pj_, 0), 



were used to rule out the enhancement of Fg from a non-uniform inclusive spectrum p{x) of 
the flnal hadrons. The obtained results have been collected in Fig. la in double logarithmic 
F2 vs M plots for different ^/s (or Mx , which is the invariant mass of the 7*Cq system 
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i.e.§ 




(4) 



in a corresponding diffractive lepton-nucleon scattering event). It is clear that the higher 
the invariant mass of 7*Cq is, the larger the power of FM's, i.e. the stronger the dynamic 
fluctuation is. In the very low invariant mass {^/s = 4.5GeV, say), the powers of FM's 
become less than 0, which can be referred to the constraint of the momentum conservation 
in the hadronization process^. 

(B). Based on pomeron-exchanged modeli, several Monte Carlo generators, such as 
POMPYTi, RAPGApli have been set up. These models have been used to describe quite 
well the HERA data in the global and averaged features, such as the rapidity gap distribu- 
tion, the diffraction and total cross section and so onlll. The theoretic models are usually 
confronted with incorrigible discrepancy, when the data about the locally nonstatistic fluc- 
tuation in small phase space are involved in the comparison with themi. No evidence has 
shown that this kind of locally dynamic fluctuation could be certainly referred only to the 
hadronization process but have nothing to do with the initial stage of high energy collisions. 
In this respect, it is also relevant to see whether the features of scale invariance and fractality 
in small phase space of the lepton-nucleon diffractive processes, specially their dependence 
upon the diffractive variables can be compatible with the pomeron type of model. In the fol- 
lowing we take RAPGAP as an example, in which the virtual photon (7*) interacts directly 
with a parton constituent of the pomeron either in lowest order (Fig.2a) or in 0(aem«s) 
order — via the photon-gluon fusion (Fig. 2b). In both cases a color octet remnant is left at 
low pt with respect to the pomeron and hence also to the proton. The higher order gluon 
emission is simulated with the Color Dipole Model (ARIADNE)§ and the hadronisation is 
performed using the JETSEtS. The pomeron flux factor /pp(t, Xp) and pomeron structure 
function G{f3) are taken as 

/pp(t,xp) = %^/3^-^-«, (5) 
given by Berger et al. and StrengS, and 

/3G(/3) = 6(1 - /?)^ (6) 
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when the pomeron is made of 2 'unrealistically hard' gluons as suggested by Ingelman 
and Schleini. Here the kinematic variables t, xp and (3 are defined as t = {P — P'Y, 
xp = '^^^.p^ ; and /3 = 2q-(p-P') ' where P' is the 4-momentum of the final state colorless 
remnant of nucleon. The pomeron Regge trajectory is given by ap{t) = 1 + e + a't, e ^ 0.085 
and the slope a' = 0.25 obtained from a fit to dat J. 

We generated 500,000 RAPGAP events, and divided the whole sample into 10 subsamples 
according to the invariable mass Mx of 7*Cq in the MC event. The scaling behaviors of FM's 
for different Mx intervals are shown in Fig. lb. The dependence of scaling behaviors of 
FM's upon Mx is similar with the result of JETSET in Fig. la, i.e. the powers increase for 
increasing scattering energy of 7*Cq. But it is noticeable that the intermittency index 02 for 
a given Mx interval is much larger in RAPGAP than that in JETSET. Having in mind that 
in Fig. la the colorless object is considered as a quark-antiquark pair formed by gluons and 
the Feymann graph of the diffractive process in this aspect is just same as shown in Fig. 2a, 
the difference between cases (A) and (B) is that the higher-order photon-parton interaction 
has been taken into account in RAPGAP (see Fig. 2b). It is understandable since the branch 
number of the parton cascading process in parton shower level is larger when the gamma 
gluon fusion is involved in Fig. 2b, while it is believed generally thati power-law behaviors in 
the color-string type of models are referred in large part to the randomly cascading process 
of parton energy in perturbative phase, so the fractal in the evolution processes with longer 
cascade branch would be stronger. 

In order to see the dependence of dynamic fluctuation upon the other diffractive variables 
as we argued above, we divided the Monte Carlo sample into 10 subsamples according to xb 
and Q^. The results of scaling behaviors of FM's are shown in Figs. 3a and b respectively. 
This example explicitly shows how the proposed method can be used to analyze the feature of 
the evolution precess of the colorless object in diffractive lepton-nucleon scattering processes. 
The significant x^-dependence of the intermittency index observed in this example not 
only shows that the hadronic final states depend very much on the average lifetime of the 
exchanged colorless object Cq if Cq can be indeed considered as a pomeron as simulated in 
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RAPGAP, but also the following: The longer the average lifetime of Cg is, the stronger the 
fractal structure of the space-time evolution processes originating from the 7*Cq's interaction 
will be. The smaller the becomes, i.e. the larger transverse dimension of virtual 7* is, 
the more 7* "behaves like a hadron". 

Since may be interpreted as the fraction of the 4-momentum of the proton carried by 
the exchanged colorless objects, and /? as the fraction of the 4-momenta of the exchanged 
objects carried by the parton interacting with the virtual boson, it meaningful to calculate 
the dependence of scaling behavior upon xp and /? (Fig. 3c and d respectively). It seems 
from RAPGAP that only when the momentum of colorless object is large enough, can scaling 
behavior of the final state originating from 7*09 be significant; and it is clear from Fig. 3d 
that the dynamic fluctuation in the diffractive scattering don't increase monotonously for 
increasing the fraction of energy of stricken parton in the pomeron. 

Last but not least, the following should be mentioned. Having in mind that jets have 
been observed (see e.g. the second paper in Ref.S and the papers cited there.) in diffractive 
electron-proton scattering processes, and HERA calorimeters have been used to measure 
transverse energies distribution of the collision events to study jet structure, it is meaningful 
to measure the scaling behavior of evolution process of 7*Cg by using the distribution of 
transverse-energies in phase space, instead of conventional multiplicity analysis. Here, the 
transverse-energy is measured on an event-by-event bases with respect to the axis of 
the virtual photon. According to the most recent experimental knowledgeB we expect to 
see that the distributions of in phase space in such collision events are symmetric with 
respect to this axis, and symmetric with respect to the origion of the c.m.s. frame of the 
colliding objects 7*Cg. As is knownll, the factorial moments defined as Eq.( of multiplicity 
of final state particles can rule out the statistical fluctuation around probability "Pm by which 
a particles appear in the m'th bin of phase space, i.e. = Cq = (1/M) Y.m=i{Pm) / {PmY- 
In order to measure the scaling behavior of probability moments of transverse energies, a 
straightforward manner following the usual procedure is to introduce an energy-unit e and 
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write the 'transverse energy factorial moment' F^^^ as {E^{E^-e) ■ ■ ■ [E^-{q-l)e]) / {E^_)i. 
It is clear that Ej_/e can be considered as integers, provided that e is sufficiently small. Under 
this condition, the statistical fluctuations around transverse energies can be removed in F^^^ 
in the same way as that in Eg defined in Eq.( H). But, this means, there is a dependence 
on an arbitrary parameter e, when we use Ejj^^l In order to check the possibility of getting 
rid of this kind of arbitrariness in the practice, let's introduce a variable, A = £/E^j_, i.e. the 
ratio between the arbitrarily chosen energy- unit e and total transverse E\_ of an event. We 
generate the transverse energies of the 'events' in the phase space by computer according to 
the Bernoulli distribution of A. It is obvious that the slope in the double logarithmic E^^'' vs 
M plot has to be flat, since there is no dynamical fluctuation in this sample. In this sample, 
we calculated the transverse energy moment i?^^-* = {E'[) / {E±)'' for different choices of A. 

( 

The corresponding R2 vs M plots is shown in Fig. 2. Here we see that the transverse energy 

( F) ( F) 

moment R2 can be considered as a good approximation for E2 , i.e. probability moments 
of transverse energies, when e, which depends upon the resolving power of the calorimeters, 
is of the order of 10"'^ of the total E±_ in the events under consideration. 

In conclusion, we have shown in this note that the colorless object Cg in diffractive lepton- 
nucleon scattering can be probed in a model-independent way, i.e. by performing the scaled 
factorial moment analyses for the hadronic flnal states originating from the scattering of 
virtual boson 7* and the Cq on an event-by-event bases. The JETSET and RAPGAP has 
been discussed as illustrative examples. It is shown in particular that higher-order inter- 
actions in hard parton level are important in the intermittency analysis of the diffractive 
scattering. Furthermore, we also pointed out that HERA calorimeter environment is a pos- 
sible place to carry out transverse-energy moment analyses of intermittency for probing 
dynamical fluctuations. 

Most of the ideas discussed in this paper were generated in conversations with Meng 
Ta-chung to whom I am grateful for patience and understanding. Thanks are also due to 
C. Boros, D. H. E. Gross, Z. Liang, R. Rittel, K. D. Schotte and K. Tabelow for helpful 
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discussions, H. Jung for the patient helps in RAPGAP Monte Carlo generator. I also thank 
Alexander von Humboldt Stiftung for the financial support. 
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FIGURES 

Fig. 1. The scaled factorial moments F2 versus the number M of subintervals of 3-dimensional 
4') phase space in log-log plot, and the second-order intermittency index 02 in corresponding 
sample set. (a). The MC result of JETSET 7.4 in different cms energy ^/s. 50000 events are 
generated in each sample set; (b). That of RAPGAP in corresponding invariable mass Mx interval 
with A'event cvents in each subsample. The pomeron flux and pomeron structure function has given 
by Eqs. ( |5|) and ( ^) in the text as a default of the MC generator. 

Fig. 2. The basic processes included in the rapgapH implementation for inelastic lepton 
scattering on a pomeron: (a) the lowest order process for hard parton level; (b) the 0{aemOis) 
order process for gamma gluon fusion. 

Fig. 3. The dependence of 2-order intermittency index (1)2 in RAP GAPS Monte Carlo imple- 
mentation upon diffractive variables: (a) x_b, (b) Q^, (c) xp and (d) p. 

Fig. 4. The 2-order transverse energy moment -^2^'' = {E'j_) / (E^)"^ as functions of partition 
number of phase space M = A/ 6 in log- log plot, when the transverse energy E± in units of e in 
subinterval 5 is stochastically produced according to Bernoulli distribution. Here, A = e/E*_^, and 
E\ is the total transverse energy in the considered phase space A. 
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